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The /3-decays of very neutron rich nuclides in the Co-Zn region were studied experimentally at the 
National Superconducting Cyclotron Laboratory using the NSCL /3-counting station in conjunction 
with the neutron detector NERO. We measured the branchings for /3-delayed neutron emission (P„ 
values) for 74 Co (18±15%), and 75_77 Ni (10±2.8%, 14±3.6%, and 30±24%, respectively) for the 
first time, and remeasured the P n values of 77_79 Cu, 79 ' 81 Zn, and 82 Ga. For 77_79 Cu and for 81 Zn 
we obtain significantly larger P n values compared to previous work. While the new half-lives for 
the Ni isotopes from this experiment had been reported before, we present here in addition the 
first half-life measurements of 75 Co (30±11 ms) and 80 Cu (170^50° ms )- Our results are compared 
with theoretical predictions, and their impact on various types of models for the astrophysical rapid 
neutron capture process (r-process) is explored. We find that with our new data the classical r- 
process model is better able to reproduce the A — 78 — 80 abundance pattern inferred from the 
solar abundances. The new data also influence r-process models based on the neutrino driven high 
entropy winds in core collapse supernovae. 



INTRODUCTION 

The rapid neutron capture process (r-process) is tra- 
ditionally believed to produce roughly half of the heavy 
elements beyond the iron region [THS]. The site of the 
r-process is still not known with certainty. The neces- 
sary very high densities of free neutrons require extreme 
conditions that have been proposed to be encountered in 
various sites within core collapse supernovae, for exam- 
ple the neutrino driven wind in delayed explosion models 
0HS], jets [7], fallback material [5], or, maybe, prompt ex- 
plosions [H1[TD]. Alternatively, the r-process could also oc- 
cur in neutron star mergers [TT] , 7-ray bursts (TJH [13] , or 
quark novae |14j . What is known from observations is the 
pattern of isotopic and elemental abundances that the r- 
process produces. This r-process abundance pattern can 
be extracted from the solar system abundances by sub- 



tracting the contributions from the s- and p-processes. 
It can also be observed directly in a specific class of ex- 
tremely metal poor but r-process element enhanced stars 
in the halo of the Galaxy (for a recent review see |15|). 
More than two dozen of such stars have been found so far 
revealing a consistent r-process pattern for elements from 
Ba to the Pt peak, but showing variations for lighter and 
heavier elements. Ongoing large scale surveys of Galac- 
tic halo stars together with high resolution spectroscopic 
followups are expected to find many more such stars in 
the future. This continuously increasing body of observa- 
tional information needs to be compared and interpreted 
with r-process models for the various sites that have been 
proposed. This requires a solid understanding of the un- 
derlying nuclear physics, that can have as much influence 
on the r-process abundances as the astrophysical environ- 
ment [T5lR8j . Currently, nuclear physics uncertainties 
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prevent the reliable extraction of site specific signatures 
from observational data and reliable calculation of the 
nucleosynthesis products of a specific r-process scenario. 

Most models assume that the r-process is a sequence 
of rapid neutron captures and /3-decays. Neutrino inter- 
actions and fission processes might play some role de- 
pending on the specific r-process environment. Among 
the most important nuclear physics quantities needed 
in r-process models are the /3-decay half-lives of the r- 
process waiting points, which determine directly the pro- 
cess timescale and the produced abundance level at their 
location in the path. Branchings for /3-delayed neutron 
emission are also important |16j during and after the r- 
process freezes out and the unstable nuclei along the r- 
process path decay towards stability. /3-delayed neutron 
emissions during that stage modify the final abundances 
and increase the neutron abundance during freezeout. In 
recent years great progress has been made in experimen- 
tally determining /3-decay properties of nuclei relevant 
for the r-process [19-28! . Nevertheless, only a very small 
fraction of the r-process isotopes have been reached ex- 
perimentally so far, most of them located near or in 
between the N = 50 and N = 82 shell closures. It- 
process model calculations therefore rely on global the- 
oretical models for the prediction of /3-decay properties 
far from stability. These models have to be tested sys- 
tematically by comparison with experimental data along 
isotopic chains far from stability. 

This collaboration has reported half-life measurements 
of the waiting point nucleus 78 Ni and other neutron rich 
Ni isotopes performed at the National Superconducting 
Cyclotron Laboratory (NSCL) at Michigan State Uni- 
versity [25]. With this measurement, the half-lives of all 
the relevant N = 50 waiting points in the r-process, 78 Ni, 
79 Cu, and 80 Zn, are known experimentally. This mass re- 
gion plays a critical role in the subset of r-process models 
(for example [10l [TBI |29] ) that are characterized by a neu- 
tron capture flow through N = 50, where it represents 
the first major bottle-neck for the production of heavier 
nuclei after the A = 8 stability gap. 

In this paper we present /3-decay half-lives and branch- 
ings for /3-delayed neutron emission for a range of very 
neutron rich Co, Ni, Cu, and Zn isotopes obtained in 
the same experiment. These data provide tests for the- 
oretical models used to predict /3-decay properties for r- 
process model calculations. In particular, the combined 
measurement of the half-life and the branching for /3- 
delayed neutron emission provides a stringent test prob- 
ing the /3-strength function at low excitation energies 
and just above the neutron threshold [3D]. In addition, 
improved data on the branchings for /3-delayed neutron 
emission of 78 Cu and 79 Cu are direct input in r-process 
model calculations. They determine the final abundance 
pattern produced by the decay of the abundances accu- 
mulated during the r-process at the N — 50 waiting point 
nuclei. 



The region around Ni is also of considerable inter- 
est for nuclear physics. In general, doubly magic nu- 
clei provide a testing ground for single particle structure 
and shell models, which in the cases of 78 Ni and 132 Sn 
are uniquely located at extreme neutron excess. Conse- 
quently, a great deal of experimental [3TH38] and theoret- 
ical [38H4T] activity has been devoted to nuclei near 78 Ni. 
An example is the existence of 8 + seniority isomers that 
is now established in 70 Ni [3TJ, 76 Ni [13137], and 78 Zn 
[4U] and is interpreted as evidence of the persistence of 
strong shell gaps out to 78 Ni. The non-existence of these 
isomers in 72 Ni and 74 Ni [34] is now explained with a 
subtle change in the residual nucleon nucleon interaction. 
Another example is the shift of the /s/2 proton orbital 
observed in the odd A Cu isotopes with increased filling 
of the neutron g 9 / 2 orbital due to the monopole term in 
the proton-neutron residual interaction [331 SI] • An un- 
derstanding of these nuclear structure effects around 78 Ni 
requires reliable shell model calculations 38, 40, 41, 44! . 
Recently some efforts have been made to develop a new 
effective shell model interaction for this region (JJ4A) 
[45] . Our new data provide a test for how well these shell 
model calculations extrapolate to the most neutron rich 
nuclei. 



SETUP 

Extremely neutron rich nuclei were produced at the 
NSCL by fragmentation of a 15 pnA 140 MeV/nucleon 
86 Kr 34+ primary beam on a 376 mg/cm 2 Be target. Af- 
ter separation with the A1900 fragment separator using 
the Bp — AE — Bp method [48], the mixed beam was 
implanted continuously into the NSCL Beta Counting 
System [45] consisting of a stack of Si PIN detectors, a 
double-sided Si strip detector (DSSD) for implantation 
of the ions, and a set of six single sided Si strip detec- 
tors (SSSD) followed by 2 additional Si PIN diodes (see 
Fig.§. 

Each implanted ion was identified event-by-event 
through a measurement of its magnetic rigidity, which 
in connection with a time-of-flight measurement provides 
the mass to charge ratio, and energy loss in the Si PIN 
detectors, which provides information about the atomic 
number. The magnetic rigidity was determined by a posi- 
tion measurement at the intermediate dispersive image of 
the A1900 fragment separator using a position sensitive 
plastic scintillator. Time-of-flight was measured between 
two plastic scintillators located at the exit of the A1900 
and in front of the Beta Counting Station (BCS). 

The resulting particle identification spectrum can be 
found in reference [25] . The implantation detector was 
a 985 pm thick DSSD with 40 x 40 pixels registering 
time and position of ion implantations and /3 decays. To 
increase the /3-decay correlation efficiency, correlations 
between decays and preceding implantations were estab- 
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FIG. 1. Schematic arrangement of the Si detectors comprising the NSCL Beta Counting System in the configuration used for 
this experiment. The degrader was adjusted to ensure implantation within the DSSD. 



fished using a 3 x 3 pixel area centered around the im- 
plantation location. The total implantation rate into all 
1600 pixels of the BCS was always less than 0.1 ions 
per second, providing enough time between implanta- 
tions into the same 3x3 pixel area for all decays of 
interest to occur long before the implantation of the next 
ion. The SSSD detectors and the downstream Si PIN 
diodes where used to veto beam contaminants that were 
not stopped in the DSSD. 

In this paper, we also present results from NERO 
(Neutron Emission Ratio Observer), a neutron detector 
that surrounded the BCS to provide information about 
branchings for /3-delayed neutron emission [5U] (see Fig. [| 
and also [26J [28]). NERO is a neutron long counter con- 
sisting of 60 ionization counters filled with BF3 or 3 He 
gas imbedded in a block of polyethylene. The counters 
are arranged in three concentric rings around the beam 
axis with the inner ring using 3 He detectors, while the 
two outer rings use BF 3 gas counters. The center of the 
detector is a cylindrical cavity with 22.4 cm diameter 
containing the beam line and the BCS, with the center 
of the DSSD being located in the center of the NERO 
detector. Neutrons emitted in the DSSD during /3-decay 
of a neutron rich isotope are thermalized in the polyethy- 
lene and then detected either by a 10 B(n,a) reaction in 
a BF 3 gas counter or a 3 He(n,p) reaction in a 3 He gas 
counter. Neutron detection events are registered using a 
multi hit TDC for each of the 60 detectors that is started 
with the detection of the /3-decay event in the DSSD. The 
typical thermalization timescale for neutrons with initial 
energies of a few MeV (for example for a 252 Cf source) 
is of the order of 70 fj,s. To determine whether a neu- 
tron has been emitted after a /3-decay we therefore count 
neutrons for a 200 fj,s time period. This is sufficient to 
detect most neutrons (about 96%) [50]. With our low 
overall implantation rate of less than 0.1 ions per second 
this is still fast enough to correlate the neutron detec- 
tion uniquely with a beta-decay event and the preceding 
implantation event. 



RESULTS FOR /3-DECAY HALF-LIVES 

The mixed rare isotope beam used in this experiment 
contained 73 " 75 Co, 75 " 78 Ni, 77 - 80 Cu, 79 ~ 81 Zn as well 
as smaller amounts of more stable isotopes 81_82 Ga, 
76 Cu, and 74 Ni. /3-decay half-lives were determined as 
described in |25j using a maximum likelihood analysis 
method (MLH) [2S]. The likelihood function to be max- 
imized is the product of the probability densities for 
each implantation event describing the measured time 
sequence of decay-type events following the implantation 
within a correlation time r corr . The formalism takes into 
account three decay generations and a constant back- 
ground. For this analysis we used r corr = 5 s. We also 
take into account /3-delayed neutron emission for the par- 
ent and daughter nuclei, which changes the daughter and 
grand-daughter half-lives. The only free parameter is the 
assumed parent half-life. Daughter half-lives and P n val- 
ues are fixed and were taken from literature [5T] and j^T] 
or, when available, from this experiment. The advantage 
of this method is that it avoids the loss of information 
about the exact time, pixel location, and time sequence 
of the multiple decay events following individual implan- 
tations of parent nuclei that occurs when binning events 
to form decay curves. 

To determine the detection efficiency for /3 particles we 
performed a traditional decay curve analysis for the cases 
with high statistics that included 75 Ni, 76 Ni, 77 Cu, 78 Cu, 
79 Zn, and 80 Zn. The half- lives obtained from the fits 
agree with the ones from the MLH analysis. The /3-decay 
detection efficiencies for the DSSD can be determined 
from the parent decay component and the known number 
of implantation events. The resulting efficiencies varied 
for different isotopes from 40.7% to 43.0%, in some cases 
by more than the statistical one sigma errors, but without 
any obvious trend. For the isotopes with less statistics 
we therefore assumed an average (3 detection efficiency of 
42±1 % for all three generations. 

The background rate for decay-type events was deter- 
mined for each implantation pixel and for each run (typ- 
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FIG. 2. Schematic arrangement of the Beta Counting Station (BCS) within the NERO neutron detector. 



ically an hour long) by counting all decay events outside 
of a 100 s time window following an implantation that 
occurred in the corresponding correlation area. Typi- 
cal background rates during the experiment ranged from 
0.008 s" 1 to 0.015 s~ 2 averaged over the entire detec- 
tor. The largest background rate was found near the 
center of the DSSD, where implantation rates are high- 
est, with a smooth decrease by about a factor of 10 to- 
wards the detector edges. This position and run depen- 
dent background rate was compared with the background 
rate determined from decay curve fits by averaging over 
the particular implantation pattern (pixel and run) for 
the individual isotope and agreed reasonably well. 

The resulting /3-decay half-lives together with the num- 
ber of implanted isotopes are listed in Table [TJ For the 
convenience of the reader we include our previously pub- 
lished half- lives for the Ni isotopes [25] . 

The errors include statistical and systematic errors. 
The statistical error is obtained from the maximum like- 
lihood analysis. To determine the systematic error we re- 
calculated half-lives for all possible variations in the input 
parameters such as efficiency, P n values, and half-lives of 
later decay generations. The envelope of all half-lives 
with their statistical errors was then used as an estimate 
for the total error. Systematic errors typically represent 
a small fraction of the total error. Contributions of the 
order of 20-30% of the total error budget are obtained for 
75 ~ 77 Ni and 77_79 Cu with parent P n values and daughter 
half-lives being the dominant sources. 



RESULTS FOR P n VALUES 

P n values were determined from the number of /3-n 
coincidences Np_ n detected within a correlation time r 
after implantation using 



p _ Np n - Njinb — Nfind /-.n 

ep and e n are j3 and neutron detection efficiencies re- 
spectively. The expected number of detected /3-n back- 
ground events Np n ^ = rg n b-^imp T corr can be calculated 
from the background rate rp^, the correlation time r corr , 
and the number of implanted ions iV; mp . The number of 
actual parent /3-decays Np = A r i, np (l — exp~ Arcorr ) can be 
determined from the known parent /3-decay rate A. Np n d 
are neutrons from /3-delayed neutron emission of daugh- 
ter nuclei, which can be significant in some of the cases 
studied here. It can be determined from 

Npnd = N imp e n ep ((1 - P n )n + P„m) (2) 

with no and n\ given by 



Pnx^K ( 1 - exp- Ar ™" 1 - exp- A - 




Indices x — 0, 1 indicate the daughters reached without 
and with neutron emission of the parent, respectively. A 
is the parent decay rate, A^ and P nx are daughter decay 
rates and P n values, respectively. Note that Np nd does 
depend on the unknown parent P n value. Eq. [T] therefore 
needs to be rearranged to calculate P n (see also f2"g]). 

The neutron detection efficiency e„ has been measured 
for a number of neutron energies in a separate experi- 
ment at the University of Notre Dame using the 13 C(a,n), 
11 B(a,n) and 51 V(p,n) reactions. In the case of 13 C(a,n) 
and 10 B(a,n) the neutron production was inferred from 
the well known properties of narrow resonances and a 
measurement of the beam current using a Faraday cup 
with electron suppression, while in the case of 51 V(p,n) 
the neutron yield was determined by offline counting of 
the induced 51 Cr activity. In addition, the detection 
efficiency was also determined with a calibrated 252 Cf 
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TABLE I. Number of detected implanted ions JV; mp , half-lives 7i/ 2 from this and previous work, number of detected correlated 
/3-n coincidences Np n , number of expected background (3 n coincidences from random background and daughter decays Np n b+d 
and P n values from this and previous work [211 [53] . 
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source. Calibration measurements using the same 252 Cf 
source were also performed before and after the exper- 
iment reported here. The details of the efficiency cali- 
bration are discussed in [50], see also [25]. The NERO 
efficiency as a function of neutron energy is constant fol- 
low energies and drops for increasing neutron energies 
beyond « 1 MeV. The energy of each neutron detection 
event is recorded to monitor correct operation of the de- 
tector system and setting of the detection thresholds, but 
does not contain information about the initial neutron 
energy because of the thermalization process. Therefore, 
theoretical assumptions about the maximum neutron en- 
ergy have to be made. From numerous past studies it is 
well known that /3-delayed neutron spectra for medium to 
heavy nuclei are compressed at energies below ~800 keV, 
often below 500 keV (for example [53], see also more de- 
tailed discussion and references in [25]). The reason is 
the phase space for /3-decay that strongly favors neutron 
emitting states just above the neutron threshold, and a 
tendency for neutron decay into excited states in the final 
nucleus. In addition, we used our shell model calculations 
[47j to estimate maximum neutron energies for 74_76 Ni, 
and 78_79 Cu. These calculations predict neutron ener- 
gies in the range of 0.9 - 1.2 MeV when assuming the 
neutron decay of the daughter states populated by /3 de- 
cay proceeds to the ground state. This assumption obvi- 
ously leads to an overestimation of the neutron energies. 
We adopt a conservative upper limit of 1.2 MeV for the 
neutron energies in this experiment, and adopt the cor- 
responding efficiency range of 33% to 41% as systematic 
error, resulting in e„=37%±4%. 

The /3-n neutron background rate r^-nb was deter- 
mined from the /3-n coincidence events among the same 



/3-type events that were used to determine the /3 back- 
ground. It was found that across the detector the ra- 
tio of /3-n to all /3 type events was constant. We there- 
fore determined this ratio for each run, and applied it 
to the position dependent /3-background rate used for 
the half-life analysis to obtain a position and run num- 
ber dependent /3-n background rate. The /3-n to /3 event 
ratio was on average about 4% on the first day of the ex- 
periment with significant scatter on an hourly timescale, 
and increased then to about 7% This is about 70 times 
larger than expected from random coincidences with the 
NERO singles neutron background rate (which was the 
same with and without beam) of 5 s _1 and the 200 /is 
/3-n correlation time window. A possible explanation is 
the light ion contamination in the radioactive beam that 
deposits energies similar to /3 particles in the DSSD and 
might emit neutrons when interacting with the detector 
stack. A majority of such events were discarded as they 
pass through the entire detector stack and were therefore 
readily vetoed, but the remainder might create a corre- 
lated /3-n background. Such an enhanced background 
was also observed in a previous NSCL experiment with 
a similar setup [26] . Background from random coinci- 
dences between actual correlated parent and daughter 
/3-decays and uncorrelated neutrons was also considered, 
but found to be negligible in all cases. 

The resulting P n values are listed together with the 
number of detected /3-n coincidences and the expected 
background neutrons in Tab. [I] In most cases the sta- 
tistical error of the number of detected neutron events 
dominates the error entirely. For the cases with high 
statistics, the other major contribution to the error is 
the uncertainty in the neutron detection efficiency. It 
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becomes comparable to the statistical error for 78 Cu and 
dominates the uncertainty for 77 Cu. Only in the case of 
74 Co there is another significant (more than 10%) con- 
tribution to the error. Here the uncertainty of the P n 
value of the daughter, 74 Ni, contributes as well. Uncer- 
tainties in background rate, /3-efficiency, number of im- 
planted ions, and daughter decay rates turn out to be 
negligible in all cases. 



DISCUSSION 

Figs. [3][7| show our new half-lives and P n values to- 
gether with previous measurements. We include the half- 
lives of the neutron-rich Ni isotopes that were published 
and discussed earlier [25]. We now present additional 
half-life data for neutron rich Co, Cu, Zn, and Ga iso- 
topes, including the first measurement of the half-lives of 
75 Co and 80 Cu. In most cases there is excellent agree- 
ment with previous measurements, except for the 79 Zn 
half-life, for which we obtain a significantly shorter value. 
It has been speculated in the past that the existence of 
an isomer in 79 Zn cannot be excluded |55j . As the popu- 
lation of isomers depends on the production mechanism, 
this could lead to differences in measured half-lives be- 
tween our technique and previous measurements using 
the ISOL method for isotope production. 

In addition to the half-lives, we also obtained P n values 
or upper limits for most of the isotopes in the beam. For 
73 ~ 75 Co, 74 ~ 77 Ni these are the first direct measurements 
of /3-delayed neutron emission. In the case of 77 ~ 79 Cu 
our measurements are more precise and systematically 
larger than previous work [55]. As stated by the authors, 
those experiments at CERN ISOLDE were difficult due 
to the chemical non-selectivity of the plasma ion source, 
and have large uncertainties. Our P n value for 79 Zn and 
81 Ga as well as the upper limits for 80 Zn and 81 Ga are 
compatible with previous work. In the case of 81 Zn, our 
P n value is larger than the literature value [56 . This is in 
line with new data from ISOLDE indicating a P n value 
of larger than 10% for 81 Zn [52]. 

Recently, Mazzocchi et al. |53j reported preliminary 
lower limits for the P D values of 73 ' 74 Co of 9±4% and 
26±9% respectively from a 7-ray spectroscopy experi- 
ment. Within their uncertainties these limits are com- 
patible with our data, albeit only barely in the case of 
73 Co, which would have to have P n between 5-8% to be 
compatible with both, our upper limit and their lower 
limit. 

After completion of this work, new experimental P n 
values for 76-78 Cu obtained by /3-delayed 7-ray spec- 
troscopy have been reported [57] (see Fig. |3j). The re- 
ported value for 76 Cu of 7.0±0.6 % is consistent with 
our upper limit. For 77 Cu a rather precise value of 
30.0±2.7 % was obtained, in excellent agreement with 
our result. For 78 Cu the reported result of 65±8 % is 



significantly larger than our measurement. This might 
reflect the difficulty of obtaining a reliable P n value via 
7-ray spectroscopy in this particular case as indicated by 
the authors [57]. 
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FIG. 3. Half-lives and P n values for the Cu isotopes mea- 
sured in this work (solid black circles) compared to previous 
work (open black triangles), and theoretical predictions from 
QRPA97 [58] (solid red line), QRPA03 g2] (dotted red line), 
CQRPA [45] (short dashed green line) and OXBASH shell 
model (see text) (long dashed blue line). The experimental 
results from [57] are added as solid black triangles. 
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In Figs. [3] - [7] our experimental results are compared 
with various theoretical calculations. The shell model 
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calculations have been performed using the JJ4B effec- 
tive interaction [22 for the P3/2/5/2P1/259/2 model space 
on top of a 56 Ni inert core. JJ4B has been developed 
starting from a realistic G-matrix interaction based on 
the Bonn-C NN potential together with core-polarization 
corrections [33]. It has been fitted to reproduce new ex- 
perimental data separately for 57_78 Ni isotopes (purely 
neutron interaction) and for N = 50 ( 79 Cu- 100 Sn) iso- 
tones (purely proton interaction). This interaction has 
been extended further to incorporate the proton-neutron 
degree of freedom. The new version of the effective in- 
teraction JJ4B [55] takes into account information about 



450 states in 73 nuclei including recent experimental data 
for Cu, Zn, Ga, Ge, As, Se, Br and Kr in the vicinity of 
78 Ni. 

The number of chosen orbitals in these calculations is 
sufficient to achieve a satisfactory description of binding 
energies, neutron separation energies and known excita- 
tion spectra for neutron-rich nuclei in the considered re- 
gion. However, we find that to reproduce the half-lives 
and P„-values of neutron-rich Ni, Zn and Cu isotopes the 
GT operator has to be renormalized by a factor of 0.37, 
instead of the expected factor of 0.75 for this mass re- 
gion. Such a strong renormalization of the GT operator 
indicates that shell model configurations responsible for 
considerable amounts of GT strength are not accounted 
for. Such configurations can be attributed to the ex- 
cluded f?/2 orbital, which is connected to its spin-orbit 
partner, f 5 / 2 , by a very strong GT matrix element. In- 
deed, a test calculation in an enlarged model space that 
includes the proton f 7 / 2 orbital and uses a 48 Ca core has 
been performed with a combined effective interaction |60j 
and confirms that the exclusion of the f? /2 orbital in the 
model space represents a strong limitation. This demon- 
strates the importance of /3-decay data, including P n val- 
ues, in testing shell model calculations far from stability. 

A global nuclear structure model is needed for astro- 
physical applications that is not limited to nuclei within 
a specific model space or near shell closures. We there- 
fore compare our data to theoretical results from the 
global, universal quasi-particle random-phase approxi- 
mation (QRPA) model the details of which can be found 
in Refs. 021 EH [SH [62] • This model considers Gamow- 
Teller /3-decay transitions from a parent nucleus to the 
accessible states in the daughter. The starting point is 
to calculate nuclear wave functions for the ground-state 



shape of the nucleus in a folded- Yukawa single-particle 
model. Ground state shapes are taken from the finite 
range droplet model (FRDM) [63]. The decay rates are 
obtained as matrix elements of the Gamow- Teller opera- 
tor between parent and daughter states in a QRPA with 
pairing and Gamow- Teller residual interactions. A global 
table of calculated /3-decay half-lives and /3-decay de- 
layed neutron emission rates was published in 1997 [58] . 
For that calculation /3-decay Q- values had been obtained 
from the 1989 Atomic Mass Evaluation [53] when avail- 
able, otherwise from the FRDM [S3] In our comparisons 
to experimental data we denote this theoretical data set 
QRPA97. 

Subsequently several enhancements have been made 
to this model resulting in a new dataset QRPA03 [32] : 
An empirical spreading was applied to the Gamow- Teller 
strength function, and for nuclei near magic numbers an 
exact spherical shape was assumed instead of the weakly 
deformed shapes obtained for these nuclei in the FRDM 
[63] . In addition, a first-forbidden strength distribution 
as predicted by the gross theory [65 , 66J was added. Com- 
pared to the allowed Gamow- Teller strength, which over 
a given energy range is represented by relatively few 
strong peaks, the first forbidden strength with its nu- 
merous small densely spaced peaks to a good approxima- 
tion constitutes a "smooth background". It is therefore 
a reasonable approach to calculate the GT transitions 
in a microscopic QRPA approach and the ff transitions 
in a macroscopic statistical model, in analogy with the 
macroscopic-microscopic method used for mass models. 
The experimental masses used to determine the Qp value 
were from [67] when available. 

We can now compare both, QRPA97 and QRPA03, 
with our new experimental data in the Co-Zn region. 
Overall QRPA03 agrees better with our set of measure- 
ments, with the exception of the Cu isotopes, where com- 
pared to QRPA97 discrepancies increased, especially for 
the /3-delayed neutron-emission probabilities (Figs. [3]). 

There are several possibilities for the origin of the dis- 
crepancies: The calculated half-life depends critically on 
the Qp. For example, for 78 Ni, a change of Qp of 1 MeV 
would change the half- life by roughly a factor of 2. As the 
half-life depends on a few low lying transitions, a 1 MeV 
change in the main low lying /3-strength just below 3 MeV 
(see Fig. [8]) would have a similar effect. P n values can be 
even more sensitive to the exact location of transitions 
in the strength distribution, in particular in a case like 
78 Ni, where there is significant strength right around the 
one neutron separation energy. Therefore, calculated P n 
values are extremely sensitive to small variations in the 
calculated strength distribution near S\ n , and no model 
can predict accurately the energy levels or strength at 
these relatively high energies. 

In addition, the strength distribution depends on de- 
formation and on whether one includes first forbidden 
transitions. This is illustrated in Table |Tlj Clearly, when 



the GT-only transitions are fast (the first three rows) 
then adding ff transitions has a small effect. Fast GT 
transitions correspond to low- lying GT strength. There- 
fore adding a small ff component at these energies has lit- 
tle effect. In the case of long half-lives (last four rows) the 
effect of including ff transitions is more substantial. One 
could perhaps argue that a standard single-particle level 
diagram for this region of nuclei reveals no obvious can- 
didates for such first forbidden decays. Fig. [8] indicates 
the main single-particle component of the strongest tran- 
sitions. However, we should recall that the residual pair- 
ing and Gamow- Teller interactions considerably change 
both the energy and wave-function structure from the 
simple single particle picture. Therefore a single-particle 
level diagram can only provide rough guidance to /3-decay 
properties. More realistic nuclear interactions might well 
yield a level structure that contains first forbidden decays 
to low-lying states. We also note that the effect of even 
fairly weak deformation can be significant. 

It is therefore noteworthy that the global model still 
reproduces measured P n values to within a factor 3. The 
good global agreement between calculated and measured 
P n and Ti/2 values is obtained although no model param- 
eter was varied to adjust the model to these experimental 
quantities. 

Nevertheless the global QRPA calculations do show a 
significant undcrprediction of the high P n values for the 
Cu isotopes found in this experiment and in [57] pointing 
to some issue in this model for this particular region. As 
discussed above, there are a number of possible reasons 
that could lead to such a deviation. One such possibil- 
ity is uncertainty in the nuclear masses. Since QRPA03, 
the masses of the Zn daughter isotopes out to 81 Zn have 
been determined with high precision Penning Trap mea- 
surements [7TJ1 IZB- As we do not mix theoretical and 
experimental masses to determine Q-values this does not 
affect the /3-decay Q-values, but it allows for a precise de- 
termination of the neutron separation energies of the Zn 
isotopes, which are needed to determine the Cu P n val- 
ues. A recalculation with these new masses lead only to 
10-30% changes compared to QRPA03, by far too small 
to explain the observed discrepancy. Furthermore, in the 
case of 76 Cu both, /3-decay Q- value and daughter neu- 
tron separation energies are known experimentally with 
keV precision in the QRPA03 calculation. Clearly, the 
underprediction of P n (with the new value of 7.3 ± 0.6 
% [57]) and half-life is already present in 76 Cu. There- 
fore, masses are unlikely to be the explanation for this 
problem. The extensive 7-ray data from [57] for 77 Cu, 
once published, might help to test calculated strength 
functions in more detail to determine the cause of this 
problem. 

Our recalculation of the QRPA predictions with up- 
dated Zn masses had significant impact on the predicted 
Zn half-lives, where the /3-decay Q-values are now known 
experimentally. For example, for 79 Zn the /3-decay Q- 
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TABLE II. Calculated /3-decay half-lives and delayed-neutron-emission probabilities for four assumptions: (1) GT transitions 
and calculated ground-state shape, (2) GT transitions and spherical ground-state shape, (3) GT+ff transitions and calculated 
ground-state shape, and (4) GT+ff transitions and spherical ground-state shape. 
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0.092 


27 


33 


25 


30 


5.34 


2.82 


6.94 


4.59 


27 


74 


0.117 


18 


26 


16 


22 


9.30 


4.75 


10.19 


7.76 


27 


75 


0.092 


15 


20 


13 


18 


7.17 


5.68 


10.41 


11.98 


28 


75 


0.058 


1053 


818 


539 


460 


9.93 


13.27 


5.86 


6.06 


28 


76 


0.050 


1064 


674 


585 


433 


25.94 


36.18 


16.97 


20.49 


28 


77 


0.050 


428 


370 


226 


207 


33.58 


35.16 


22.80 


25.17 


28 


78 


0.025 


371 


371 


228 


224 


38.26 


59.53 


28.69 


39.06 


29 


76 


0.117 


734 


463 


435 


318 


8.06 


1.28 


4.74 


2.46 


29 


77 


0.083 


827 


417 


505 


318 


32.81 


10.40 


19.30 


10.38 


29 


78 


0.075 


386 


289 


224 


187 


34.82 


11.16 


19.79 


10.60 


29 


79 


0.050 


391 


224 


223 


155 


46.32 


18.97 


38.07 


25.24 


29 


80 


0.075 


103 


160 


64 


84 


58.46 


100.00 


51.53 


56.62 


30 


79 


0.067 


4991 


3591 


1869 


1647 


1.29 


1.76 


0.29 


0.34 


30 


80 


0.042 


3796 


2505 


1580 


1259 


20.91 


4.04 


5.72 


6.19 


30 


81 


0.075 


707 3160 


325 


517 


19.53 


43.42 


11.34 


13.38 


31 


82 


0.083 
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value increases from 8.68 MeV as predicted by the FRDM 
to 9.08±0.1 MeV, decreasing the predicted half-life from 
1.6 s to 1.0 s. Similarly, the 80,81 Zn half-lives decrease 
from 1.3 s, and 0.52 s to 1.0 s and 0.35 s, respectively. 
Overall this reduces the discrepancy between experiment 
and theoretical prediction significantly. 

Finally, we compare our new data with the contin- 
uum QRPA (CQRPA) calculations from Borzov 05]. His 
treatment is limited to spherical nuclei and is not global, 
but does include a pn interaction in the particle-particle 



channel and a microscopic calculation of the first for- 
bidden strength. Table |H| shows that the assumption of 
spherical shapes may be inappropriate as even the mod- 
est deformations predicted for his region are expected 
to have a substantial effect on half-lives and P n values. 
The predictions of half-lives and P n values by Borzov 
show overall very good agreement with our experimental 
results. Borzov finds, that in his model first forbidden 
transitions play only a minor role around 78 Ni as long as 
N < 50. However, the CQRPA model does significantly 
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underpredict the P n values of Cu and Cu, while there 
is excellent agreement for 78 Cu and 79 Cu. 

Motivated by their new experimental P n values for 
77-78q u Winger e t a j_ [57] ar g Ue that this underpre- 
diction is an indication for the inversion of the 7r2p 3 / 2 
and 7rl/ 5 /2 single particle orbitals. CQRPA model cal- 
culations where these orbitals are inverted indeed lead to 
larger P n values, though no half-life data are presented 
that would allow to verify consistency with experimen- 
tal half-lives. However, our P D value for 78 Cu is signif- 
icantly lower than the experimental value reported by 
|57j . While the old CQRPA results were in agreement 
with our measurement, the new CQRPA value of 53% 
is somewhat high. In addition, the old CQRPA results 
are also in excellent agreement with our new P n value 
for 79 Cu. Nevertheless, there is other experimental evi- 
dence for such a level inversion to occur for Cu isotopes 
at 75 Cu and beyond [75] thereby justifying the modifica- 
tions of the CQRPA model. 

R-PROCESS CALCULATIONS 

With our new data and recent precision mass measure- 
ments around 80 Zn [70l [71] the nuclear physics needed to 
model the r-process around A = 80 is now to a large ex- 
tent experimentally determined. We can therefore test 
r-process models in this particular mass region against 
observations with greatly reduced nuclear physics uncer- 
tainties. 

It is quite challenging to understand the origin of the 
elements in this mass region, as not only all major neu- 
tron capture processes, the weak and strong s-process 
and the r-process, can contribute, but charged particle 
reaction sequences can reach this mass region as well. In- 
deed, one class of r-process models, the neutrino driven 
wind scenario in core collapse supernovae, predict that 
nuclei in this region are produced by a combination of 
charged particle and neutron induced reactions. Never- 
theless, we can ask, whether r-process models character- 
ized by a neutron capture and /3-decay reaction sequence 
in the A = 80 mass region are now able to reproduce 
the observed solar r-process abundances in this region, 
or not. 

To address this question we use a classical r-process 
model that simulates a series of neutron exposures of Fe 
seed nuclei with neutron density n„, temperature T, du- 
ration r, and weight ui. It has been shown that such a 
model can reproduce the observed solar abundance pat- 
tern reasonably well employing power law relationships 
w(n ra ) = a\n°£ and r(n„) = a^n^ leaving only 3 free pa- 
rameters plus an overall normalization |73j . The temper- 
ature T = 1.35 GK is kept constant and is the same for 
all components representing a typical r-process freezeout 
temperature. The model we use adopts the waiting point 
approximation and assumes a sudden freezeout with de- 



cay back to stability once all exposures have been applied. 
The fact that the site of the r-process is still unknown and 
that a wide range of scenarios have been proposed, moti- 
vates the use of this simple site-independent model. The 
classical model can simulate well the final local neutron 
capture flow in an r-process scenario at freezeout, which 
tends to dominate the features of the final abundance 
pattern. Of course further modifications of the abun- 
dance pattern can occur during freezeout, but this effect 
is highly model dependent. It therefore makes sense to 
explore the agreement of a simple r-process model with 
observations. Once the nuclear physics is fixed, major 
disagreements with observations might indicate an en- 
tirely different r-process mechanism for the mass region 
in question, while smaller deviations might reveal addi- 
tional, site specific effects, such as an extended freezeout. 
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FIG. 9. Abundances calculated with a classical r-process 
model using previously known decay data (dotted black line) 
and the new data from this work (thick solid red line). The 
thin solid blue line indicates a calculation with the large r8 Ni 
half-life predicted by QRPA97. The thick dotted red line 
shows results using the new half-lives from our work, but the 
old P n values. Solar r-process abundances are shown as black 
data points. 
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Fig. [9] shows results from an r-process calculation, 
where the four model parameters have been htted to re- 
produce the solar r-process residuals from Arlandini et 
al. [71] . These r-process residuals have been obtained by 
subtracting from the observed solar abundances a main 
s-process component calculated with a classical s-process 
model and a weak s-process component. In the A = 80 
mass region this solar r-process abundance distribution is 
very similar to the one obtained with a set of realistic stel- 
lar s-process models and updated solar abundances [75] , 
Our baseline model calculation with model parameters 
ai = 4224.74, a 2 = -0.082, a 3 = 788.86 and a 4 = 0.0089 
uses the nuclear masses from the 2003 Atomic Mass Eval- 
uation |76j . including their extrapolations, updated with 
the new masses in the Zn region [TDJ, [71] . Masses of more 
exotic nuclei are taken from the ETFSI-Q mass model, 
which has been shown to be well suited for r-process 
calculations. We also performed calculations using the 
FRDM mass model that lead to very similar results in 
the A = 80 mass region and to the same conclusions. Ex- 
perimental /3-decay data, including P n values, are taken 
from NNDC [77] when available, otherwise our theoreti- 
cal QRPA03 data are used. For comparison we then run 
the same calculation with our updated /3-decay proper- 
ties from this study. We also show a calculation using 
the large 78 Ni half-life predicted in our earlier QRPA97 
model. 

The global effects of implementing our new experimen- 
tal data are rather modest. However, the half- life of 
78 Ni clearly plays an important role, and the large 78 Ni 
half-life of QRPA97 leads to noticable changes across all 
masses, including the A = 195 peak region. Our mea- 
surement of the half-life of 78 Ni now excludes such a long 
half-life with certainty. 

In the A — 80 region changes are more significant. 
Our new data lead to a more pronounced odd-even effect 
for A = 78, 79, 80 that agrees better with observations. 
As Fig. [9] shows, our shorter 78 Ni half-life actually leads 
to a weaker odd-even effect between A — 78 and 79, 
but the larger P n values for the Cu isotopes, especially 
79 Cu, more than compensate and are therefore important 
to obtain a better fit. The need for large P n values, 
especially for 79 Cu in order to reproduce the pronounced 
odd-even effect in the r-process abundances has already 
been pointed out in earlier studies based on theoretical 
predictions of half- lives and P n values [75] . 

Overall, with the masses and decay data for A < 81 
now available the A = 80 region can be reproduced quite 
well. The agreement is very good for A = 78,79,80 
where some experimental data are now available. On 
the other hand, some problems become apparent in the 
A = 81 — 90 mass range. Because of the lack of experi- 
mental constraints on the nuclear physics in this region, 
a nuclear physics explanation cannot be excluded yet. In 
this region, most of the relevant masses, and most of the 
relevant /3-decay half-lives beyond A = 84 are unknown. 



Clearly more experimental work needs to be done to ex- 
tend the mass region of reliable nuclear physics over the 
entire A = 80 abundance peak area. 

Even in the A — 78 — 80 area there are still some 
remaining nuclear physics uncertainties. The most im- 
portant nuclear physics data at the N = 50 shell closure 
are the half-lives of the major N = 50 waiting points, 
78 Ni, 79 Cu, and 80 Zn and the nuclear masses for the 
isotopic chains, where the neutron capture flow might 
begin to cross the N = 50 shell closure. Owing to the 
rather large drop in neutron separation energy across the 
N = 50 shell gap it is clear that regardless of the mass 
model adopted, for the Ni and Cu isotopes the neutron 
densities required for the reaction flow to cross N = 50 
(n n > 10 24 ' 26 for the Cu, Ni isotopic chains, respectively) 
are much higher than the neutron densities that produce 
the bulk of the A = 80 isotopes. Therefore, for the cal- 
culation of abundances around A — 80, precision masses 
to characterize the breakout beyond N — 50 in detail are 
mainly needed for Zn and Ga. For these isotopic chains 
masses are needed across N = 50 out to N — 52 because 
of the odd-even effect of the neutron separation energies. 
In addition, P n values are needed for all A = 78 — 80 
waiting points and their decay daughters. With our data 
and previous work, all of these quantities are now ex- 
perimentally known, except for the mass of 82 Zn, which 
might introduce some uncertainty concerning the 80 Zn 
waiting point. This has been discussed in detail in [70j . 
In addition, the measured 78 Ni half-life still has a large 
error bar. While the rather long half-life predicted by 
QRPA97 is now excluded, the half-life change that leads 
to the changed A = 78 abundance indicated in Fig. [9] 
is actually of the order of the experimental uncertainty. 
Therefore, a more precise 78 Ni half-life would be desir- 
able. The P n value of 78 Ni might have some influence on 
the result as well, and should be determined experimen- 
tally. 

Finally, we also explore the role of the nuclear physics 
in the 78 Ni region in a more site specific r-process model, 
the so called high entropy wind (HEW) scenario [4"H6"II79T- 
155] . The model is inspired by the conditions expected 
near the proto-neutron star forming in a core collapse 
supernova shortly after the explosion. The high neu- 
trino luminosity is thought to drive outflows of strongly 
heated, low density matter (hence high entropy) that at 
late times become neutron rich. While this is one of the 
most promising candidates for an r-process scenario, re- 
alistic models based on conditions obtained in current su- 
pernova models do not lead to a full r-process. We there- 
fore use the simplified parametrized model of [SU] [H3] 
that follows a set of one dimensional adiabatic expan- 
sions (components), each characterized by an entropy per 
baryon S, an initial electron abundance Y e , and an ex- 
pansion velocity v. Nuclear reactions are followed with 
a full reaction network including /3-decay properties and 
all neutron, proton, and 7 induced reactions, but neglect- 
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ing neutrino interactions and fission. One possible choice 
of parameters that lead to a successful r-process is to 
keep the same realistic values of Y e = 0.45 and v = 7500 
km/s for all components, but choose an equidistant set 
of entropies that ranges up to about S/k ~ 250, larger 
than predicted by supernova models. When assuming 
that equal amounts of material are processed by each 
component, such a model has been shown to reasonably 
reproduce the solar r-process residuals [501 IH3] • We use 
the same /3-decay rates and P n values employed in our 
classical model calculations. Neutron and charged parti- 
cle induced reaction rates are taken from NON-SMOKER 
statistical model predictions [84] using the FRDM mass 
model. 7 induced reactions are calculated from their in- 
verse capture reactions via detailed balance. 

In HEW models, r-process seed nuclei are produced by 
combined charged particle and neutron induced processes 
in the A ~ 90 region close to stability. As this is already 
beyond the 78 Ni region, one might expect that /3-decay 
properties near 78 Ni do not play a role. Our calculations 
show, however, that at high entropies around S/k = 200, 
where neutron to seed ratios become high enough to pro- 
duce the heavier r-process elements, 78 Ni becomes part of 
the r-process path. Fig. 10 shows the abundance distri- 



0.0001 r 



bution for S/k = 200 with the previous nuclear database, 
and with our new experimental results. It turns out that 
the change is entirely due to our new, shorter 78 Ni half- 
life (110 ms instead of 224 ms). For comparison we also 
show a calculation with the older, longer 78 Ni half-life 
(477 ms) predicted by [55]. Clearly a long 78 Ni half-life 
reduces the A — 78 production, but increases the produc- 
tion of A = 100 - 120 nuclei at S/k ~ 200. Interestingly, 
the production of very heavy r-process nuclei is slightly 
suppressed by a shorter 78 Ni half-life, contrary to what 
one would expect naively and opposite from the behavior 
in the classical model. Nevertheless, the impact of the 
78 Ni half-life on the final abundances, once all entropy 
components have been added up, is rather small. Fig.[TT] 
shows the ratio of our new abundances to what one ob- 
tains with the long 78 Ni half-life from [55]. Besides the 
significant change at A = 78, there is a 10% increase at 
A = 130 (where S/k =200 makes its largest contribution) 
and a 10% suppression of very heavy nuclei. This mod- 
est sensitivity of the final abundances reflects the rather 
narrow entropy range between S/k ~ 190 — 210 that is 
broadly influenced by 78 Ni. Effects for higher entropies 
are still significant but only below A ~ 130 where those 
high entropy components do not contribute much. How- 
ever, it is likely that critical waiting points such as 78 Ni in 
the case of S/k ~ 200 do exist also for the other entropy 
components. It will be important to identify and mea- 
sure these waiting points to obtain more reliable HEW 
r-process calculations. In addition to the already known 
N = 50 isotopes 80 Zn and 79 Cu, our first measurement 
of the 78 Ni half-life is an important further step towards 
this goal. 
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FIG. 10. Abundances calculated with the HEW r-process 
model for a single entropy component with S/k =200. Results 
based on the data from this work ( 78 Ni half-life of 110 ms) 
(thick solid red line) are compared with results based on pre- 
viously available data ( 78 Ni half-life of 224 ms) (thin solid 
blue line) and previously available data with the long 78 Ni 
half-life of 477 ms from |58j (black dotted line) . 
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FIG. 11. Ratio of final HEW abundance distributions 
summed over all entropy components calculated with our new 
data including a shorter 78 Ni half-life of 110 ms to what one 
obtains using previously available data and the long 78 Ni half- 
life of 477 ms from l58l. 



Using our new P n values instead of the old ones does 
not lead to significant changes in the calculated abun- 
dances for our choice of HEW model parameters. In 
principle one can expect a reduced impact of /3-delayed 
neutron emission compared to the classical model, as neu- 
trons are present at later times potentially reversing the 
effects via neutron capture. Once the neutrons are ex- 
hausted, the r-process path tends to be already closer 
to stability where P n values are smaller. To explore the 
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impact of P n values in HEW calculations, we run a sim- 
ulation without any /3-delayed neutron emission for com- 
parison. This is similar to what has been done in [16] 
for the classical r-process model. The result is shown in 
Fig. [12] and demonstrates the importance of /3-delayed 
neutron emission in HEW models. For A < 110 the im- 
pact is less pronounced, though for some mass chains 
significant changes of up to a factor of 2 do occur. On 
the other hand, P n values play a critical role in shaping 
the A = 130 and A — 195 abundance peaks as well as 
the rare earth peaks. 




100 150 200 

mass number 

FIG. 12. HEW abundance distributions summed over all en- 
tropy components using our new data (thick solid red line) 
and setting all P n values to zero (thin solid blue line) com- 
pared with the solar r-process residuals (black filled circles). 

Overall the impact of our new data on the synthesis 
of A ~ 80 r-process nuclei in our HEW model is rather 
modest, in contrast to what we found in the classical 
model. However, as has been discussed extensively in 
[80] . and as can be seen in Fig. 12 our particular choice of 
HEW model parameters, especially our initial Y e , clearly 
does not reproduce the solar r-process residuals below 
A ~ 110. Once this problem has been solved, the ques- 
tion of the relevance of the nuclear physics in the A = 80 
region in HEW models should be revisited. 



SUMMARY 

We have provided the first experimental P n values for 
extremely neutron rich Ni isotopes out to ' ' Ni and a first 
half-life for 80 Cu. Our experimental P n values for 77 Cu 
and 78 Cu confirm the finding of an experiment carried out 
in parallel to this study [57] that the previously reported 
P n values are too small, though in the case of 78 Cu we 
find the discrepancy is somewhat smaller as reported in 
[57j . This discrepancy is not too surprising given the 
large error bars and the difficulties encountered in these 



earlier measurements. 

Given these discrepancies with earlier work, our mea- 
surement of the P n value of 79 Cu is of particular im- 
portance because of the key role this nucleus plays in 
r-process models. We find that in this case our new mea- 
surement agrees well with work reported earlier indicat- 
ing that the problems with earlier experiments are not 
necessarily of general nature. 

Our experimental data provide constraints for theo- 
retical models used to understand the nuclear structure 
of the 78 Ni region. We find that recent shell model cal- 
culations as well as local continuum QRPA calculations 
describe the experimental data well. This gives some 
confidence that such models can be used to investigate 
nuclear structure effects around 78 Ni. However, in the 
case of the shell model calculations the required large 
GT quenching factor indicates the need for larger model 
spaces. 

Our results also show that the global QRPA03 model 
used to predict nuclear structure input for astrophysi- 
cal r-process calculations agrees with data within the ex- 
pected theoretical uncertainty. The model exhibits some 
local deficiencies in the 78 Ni mass region, in particular for 
the P n values of the neutron rich Cu isotopes. More ex- 
perimental data are needed to identify the nature of this 
problem. With our measurements experimental decay 
data are now available for r-process calculations along 
the entire reaction path at TV = 50. This includes our 
improved data on the branchings for /3-delayed neutron 
emission of 78 Cu and 79 Cu, which are needed in r-process 
model calculations to reliably calculate the final abun- 
dances in the A = 78 — 79 mass range. Together with 
recent mass measurements around 80 Zn this now puts a 
three isotope section of the r-process around A — 80 on a 
fairly solid experimental basis. Our r-process model cal- 
culations demonstrate that this narrow mass region from 
A = 78 — 80 turns out to be well reproduced with a neu- 
tron capture flow based r-process. More experimental nu- 
clear physics data beyond A = 80 are needed to broaden 
the mass range in which r-process models can be tested 
reliably, but for now neutron captures under typical r- 
process conditions cannot be excluded as a mechanism 
for the origin of the elements around A — 80 not made 
by the s-process. 

We also show that the 78 Ni half-life does play a role 
in HEW r-process models. On the other hand, compared 
to the classical model, the HEW model seems to be less 
sensitive to P n values in the Ni-Cu region, at least for 
our choice of parameters. We demonstrated that the P n 
values of heavier nuclei do play a critical role in our HEW 
r-process model. 
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